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Abstract. Properties of quantum shape-phase transitions in finite nuclei are considered in the 
framework of the interacting boson model. Special emphasis is paid to the dynamics at the critical- 
point of a general first-order phase transition. 
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Phase transitions associated with a change of shape are known to occur in dynamical 
systems such as nuclei. Recently, it has been recognized that such quantum shape-phase 
, transitions are amenable to analytic descriptions at the critical points 01 0]. For nu- 

clei these analytic benchmarks of criticality were obtained in the geometric framework 
of a Bohr Hamiltonian for macroscopic quadrupole shapes. In particular, the E(5) [Q]] 

O . (X(5) [2]) benchmark is applicable to a second- (first-) order shape-phase transition 

between spherical and deformed y-unstable (axially- symmetric) nuclei. Empirical evi- 
dence of these benchmarks have been presented y|, Hi- An important issue concerning 
phase transitions in real nuclei is the role of a finite number of nucleons. This aspect 
can be addressed in the algebraic framework of the interacting boson model (IBM) dfj] 

+f ; which describes low-lying quadrupole collective states in nuclei in terms of a system 

of N monopole (s) and quadrupole (d) bosons representing valence nucleon pairs. The 
three dynamical symmetry limits of the model: U(5), SU(3), and 0(6), describe the dy- 
namics of stable nuclear shapes: spherical, axially-deformed, and y-unstable deformed. 
A geometric visualization of the model is obtained by an intrinsic energy surface defined 
by the expectation value of the Hamiltonian in the coherent (intrinsic) state da [ZD 

\p,r>#) = (ivo-V2(&ty>) , (i) 

where b\ = ( 1 + /3 2 ) -1 / 2 [j8 cos ydj + /3 sin y (d\ + d^_ 2 ) / \fl + / ] . For the general IBM 
Hamiltonian with one- and two-body interactions, the energy surface takes the form 

E(P,y) = N(N-l)(l+p 2 )- 2 [a/3 2 -£/3 3 cos3y+c/3 4 ] . (2) 

The coefficients a, b, c involve particular linear combinations of the Hamiltonian's pa- 
rameters [[sLIsl]- The quadrupole shape parameters in the intrinsic state characterize the 
associated equilibrium shape. Phase transitions for finite N can be studied by an IBM 
Hamiltonian involving terms from different dynamical symmetry chains flg]. The nature 
of the phase transition is governed by the topology of the corresponding energy surface. 
In a second-order phase transition, the energy surface is y-independent and has a single 
minimum which changes continuously from a spherical to a deformed y-unstable phase. 
At the critical point, a = b = 0, and the energy- surface acquires a flat behaviour (~ /3 4 ) 
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FIGURE 1. The IBM energy surface, Eq. (3), at the critical point of a first-order phase transition. The 
position and height of the barrier are j3+ = and h = /(j3+) = \ 



1 + J 1 + J3q I respectively. 



for small /3. This is the situation encountered in the U(5)-0(6) phase transition. In a 
first-order phase transition the energy surface has two coexisting minima which become 
degenerate at the critical point. The U(5)-SU(3) phase transition is a special case of this 
class, for which the barrier separating the spherical and axially-deformed minima is ex- 
tremely small, and hence the critical energy surface is rather flat. Numerical |3j,LL0|] and 



(approximate) analytic 1UJ,LL2D studies within the IBM, show that the U(5)-0(6) [U(5)- 
SU(3)] critical Hamiltonians capture the essential features of the E(5) [X(5)] models 
which employ infinite square-well potentials. In the present contribution we consider 
the properties of a general first-order phase transition with an arbitrary barrier. In this 
case the critical energy surface satisfies b > and b 2 = Aac, and for y = has the form 

EcniP) = cN(N-l)f(P) 
f(P) = /3 2 (l+/3 2 )- 2 (/3 -/3 ) 2 . (3) 

As shown in Fig. 1, E cr i(fi) exhibits degenerate spherical and deformed minima, at 
/3 = and /3 = /3o = t2 > 0. The value of /3o determines the position (/3 = j8 + ) and 
height (h) of the barrier separating the two minima in a manner given in the caption. To 
construct an Hamiltonian with such an energy surface, it is advantageous to resolve the 
critical Hamiltonian into intrinsic and collective parts [a M\ 

H = H mt + H c . (4) 

The intrinsic part of the Hamiltonian (Hi„ t ) is defined to have the equilibrium condensate 
| /3 = /3o, 7= 0;N) as an exact zero-energy eigenstate and to have an energy surface as 
in Eq. (3). It can be transcribed in the form 

= h 2 P}-h, (5) 
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FIGURE 2. Left potion: spectrum of H int , Eq. (5), with h 2 = 0.1, j3o = 1.3 and N = 10. Right portion: 
the number of d bosons (n c i) probability distribution for selected eigenstates of //,„,. 



with ^ = j8o + A/5 (d^) ^ . By construction H int has a set of solvable deformed 
eigenstates with energy E = 0, which are the states |/3 = Pq;N,L) with angular momen- 
tum L = 0, 2, 4, . . . , 2N projected from the intrinsic state |/3 = /3o, y = 0; N) . It has also 
solvable spherical eigenstates: |A^,n^ = T = L = 0) = 1^) and |A^, = % = 3,L = 3) 
with energy £ = and E = 3/i2 [Pq(N — 3) +5] respectively. For large N the spectrum 
of Hi nt is harmonic, involving 5-dimensional quadrupole vibrations about the spherical 
minimum with frequency £, and both /3 and y vibrations about the deformed minimum, 
with frequencies ep and e y given by 

? 9 
£ = e p =h 2 p$N , £ y =— - ^ep. (6) 

1 +Po 

For the acceptable range < /3o < 1.4, the y-band is expected to be considerably higher 
than the /3-band. All these features are present in the exact spectrum of Hi nt shown 
in Fig. 2, which displays a zero-energy deformed (K = 0) ground band, degenerate 
with a spherical [rid = 0) ground state. The remaining states are either predominantly 
spherical, or deformed states arranged in several excited K = bands below the y band. 
The coexistence of spherical and deformed states is evident in the right portion of Fig. 2, 
which shows the decomposition of wave functions of selected eigenstates of Ht nt . The 
"deformed" states show a broad distribution typical of a deformed rotor structure. The 
"spherical" states show the characteristic dominance of single components that one 
would expect for a spherical vibrator. In particular, as mentioned above, the solvable 
L = O2 and L = 3^ states are pure U (5) states. 

The collective part (H c ) of the full Hamiltonian, Eq. (4), is composed of kinetic terms 
which do not affect the shape of the energy surface. It can be transcribed in the form 



H c = c 3 [£ ( 3 ) -6n d ] +c 5 [C 0(5) -Ah d ] +c 6 



C 0(6)- 5N 



(7) 



where h c i and N — nd + n s are the J-boson and total-boson number operators respectively. 
Here Qj denotes the quadratic Casimir operator of the group G as defined in |0]. In 
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FIGURE 3. The spectrum of the critical Hamiltonian H = Hi nt +H C , Eqs. (5) and (7), with j3o = 1.3, 
N = 10, and C3/I12 — 0.05 or C5/I12 =0.1 or c^jh^ = 0.05. The spectrum of a rigid rotor and the X(5) 
model is shown for comparison. Energies are in units of the first excited state, E(2^) —E(0^) = 100, and 
B(E2) values are in units of B(E2; 2+ -> 0j") = 100. 



general H c splits and mixes the states of Hi nt . Fig. 3 shows the effect of different 
rotational terms in H c , with parameters indicated in the caption. For the high-barrier 
case considered here, (/3o = 1 .3, h = 0. 1), the calculated spectrum resembles a rigid-rotor 
(E ~ aj^L(L+ 1)) for the C3-term, a rotor with centrifugal stretching (E ~ a^L{L-\- 1) — 
bft[L(L+ l)] 2 ) for the C5-term, and a X(5)-like spectrum for the C6-term. In all cases the 
B(E2) values are close to the rigid-rotor Alaga values. This behaviour is different from 
that encountered when the barrier is low, e.g., for the U(5)-SU(3) critical Hamiltonian 
with /3o = 1 /2y/2 and h ta 10~ 3 , where both the spectrum and E2 transitions are similar 
to the X(5) predictions. Considerable insight of the underlying structure at the critical 
point is gained by examining the 2x2 mixing matrix between the spherical (\s N )) 
and deformed (/3;A/,L = 0)) states. The derived eigenvalues serve as eigenpotentials, 
and the corresponding eigenvectors are identified with the ground- and first-excited 
L = states. The deformed states |j8;iV,L) with L > are identified with excited 
members of the ground-band with energies given by the L-projected energy surface 

E^ = (j5;N,L\Hi nt +H c \j5;N,L), which can be evaluated in closed form 
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M/3-/3o) 2 4^+c 3 \l(L+1)-6dW 



+ C 5 



(AO 
2,L 



+ C6 



ys^ + N(N-V 



(8) 



Here E 



(N) 



-2 L'^1L'^2 L an( * ^2 l denote the expectation values of h s h c \, n^, fid{nd ~ 1) an d 
n s (n s — 1) respectively in |/3;A/,L). The value of /3 in the indicated wave functions and 
energies, is chosen at the global minimum of the lowest eigenpotential. This procedure 
yields an excellent approximation to the structure of yrast states at the critical point, 
which is then used to derive accurate finite-N estimates to the corresponding energies 
and E2 rates [13]. The same prescription is applicable also to a first-order phase transi- 
tion with a low-barrier For a second-order critical point, the determination of the 
effective /3 -deformation involves an 0(5) projection without two-level mixing ifllll . 
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